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Abstract 

Full counting statistics of electron transport is a powerful diagnostic tool for probing the na¬ 
ture of quantum transport beyond what is obtainable from the average current or conductance 
measurement alone. In particular, the non-Markovian dynamics of quantum dot molecule plays 
an important role in the nonequilibrium electron tunneling processes. It is thus necessary to un¬ 
derstand the non-Markovian full counting statistics in a quantum dot molecule. Here we study 
the non-Markovian full counting statistics in two typical quantum dot molecules, namely, serially 
coupled and side-coupled double quantum dots with high quantum coherence in a certain param¬ 
eter regime. We demonstrate that the non-Markovian effect manifests itself through the quantum 
coherence of the quantum dot molecule system, and has a significant impact on the full count¬ 
ing statistics in the high quantum-coherent quantum dot molecule system, which depends on the 
coupling of the quantum dot molecule system with the source and drain electrodes. The results 
indicated that the influence of the non-Markovian effect on the full counting statistics of electron 
transport, which should be considered in a high quantum-coherent quantum dot molecule system, 
can provide a better understanding of electron transport through quantum dot molecules. 
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Introduction 


Full counting statistics P (FCS) of electron transport through mesoscopic system has 
attracted considerable attention both experimentally and theoretically because it can pro¬ 
vide a deeper insight into the nature of electron transport mechanisms, which cannot be 
obtained from the average current [2Hin]. For instance, the shot noise measurements can 
be used to probe the dynamical in an open double quantum dots (QDs) |TT], the coherent 
coupling between serially coupled QDs [12] , the evolution of the Kondo effect in a QD [13] , 
and the conduction channels of quantum conductors [Tl|. In particular, shot noise char¬ 
acteristics can provide information about the feature of the pseudospin Kondo effect in a 
laterally coupled double QDs [I5], the spin accumulations in a electron reservoir [T6|, and 
the charge fractionalization in the v = 2 quantum Hall edge DU- In addition, the degree 
of entanglement of two electrons in the double QDs [T8|, the dephasing rate in a closed QD 
[T9|, the internal level structure of single molecule magnet [201 [21] can be characterized by 
the super-Poissonian shot noise . 

On the other hand, the quantum coherence in coupled QD system, which is characterized 
by the off-diagonal elements of the reduced density matrix of the QD system within the 
framework of the density matrix theory [22]. plays an important role in the electron tunnel¬ 
ing processes and has a significant influence on electron transport [23I133] . In particular, 
theoretical studies have demonstrated that the high-order cumulants, e.g., the shot noise, 
the skewness, are more sensitive to the quantum coherence than the average current in the 
different types of QD systems [I2l IMIIHS] and the quantum coherence information in a side- 
coupled double QD system can be extracted from the high-order current cumulants [35] . In 
fact, the non-Markovian dynamics of the QD system also plays an important role in the 
non-equilibrium electron tunneling processes. However, the above studies on current noise 
or FCS were mainly based on the different types of Markovian master equations. Although 
the influence of non-Markovian effect on the long-time limit of the FCS in the QD systems 
has received some attention [33l |39]-06], how the non-Markovian effect affects the FCS is 
still an open issue, especially the influence of the interplay between the quantum coherence 
and non-Markovian effect on the long-time limit of the FCS has not yet been revealed. 

The aim of this report is thus to derive a non-Markovian FCS formalism based on the 
exact time-convolutionless (TCL) master equation and study the influences of the quantum 
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coherence and non-Markovian effect on the FCS in QD molecule systems. It is demon¬ 
strated that the non-Markovian effect manifests itself through the quantum coherence of 
the considered QD molecule system, and has a significant impact on the FCS in the high 
quantum-coherent QD molecule system, which depends on the coupling of the considered 
QD molecule system with the incident and outgoing electrodes. Consequently, it is neces¬ 
sary to consider the influence of the non-Markovian effect on the full counting statistics of 
electron transport in a high quantum-coherent single-molecule system. 

Results 

We now study the influences of the quantum coherence and non-Markovian effect on 
the FCS of electronic transport through the QD molecule system. In order to facilitate 
discussions effectively, we consider three typical QD systems, namely, single QD without 
quantum coherence, serially coupled double QDs and side-coupled double QDs with high 
quantum coherence in a certain parameter regime (see Fig. 1). In addition, we assume 
the bias voltage (pl = —Pi? = 14/2) is symmetrically entirely dropped at the QD-electrode 
tunnel junctions, which implies that the levels of the QDs are independent of the applied 
bias voltage even if the couplings are not symmetric, and choose meV as the unit of energy 
which corresponds to a typical experimental situationj47j. 

Single quantum dot without quantum coherence 

In this subsection, we consider a single QD weakly coupled to two ferromagnetic elec¬ 
trodes. The Hamiltonian of the considered system is described by the Htotal = Hdot + Hieads + 
Ht- The QD Hamiltonian Hdot is given by 

hfdotjl ad\dfj + U d\.d^d}^di, (1) 

(J 

where dj. [da] creates (annihilates) an electron with spin a and on-site energy Eo- (which can 
be tuned by a gate voltage Vg) in this QD system. U is the intradot Coulomb interaction 
between two electrons in the QD system. 

The relaxation in the two ferromagnetic electrodes is assumed to be sufficiently fast, so 
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that their electron distributions can be described by equilibrium Fermi functions. The two 
electrodes are thus modeled as non-interacting Fermi gases and the corresponding Hamilto¬ 
nians can be expressed as 

.^LeadSjl (2) 

Q:k5 

where (flaks) creates (annihilates) an electron with energy spin s and momen¬ 
tum k in a (a = L, R) electrode, and s = -|- (—) denotes the majority (minority) 
spin states with the density of states ga,s- The polarization vectors pl (left lead) and 
P/j (right lead) are parallel to each other, and their magnitudes are characterized by 
Pq = |Pa| = {9a,+ — 9a-) / {.9a,+ + 9a-)- The tuuneliug between the QD and the electrodes 
is described by 

i^T,! ~ ^Lfc+fl^fc-l-C^t T T T T H.C., (3) 

where spin-up ^ and spin-down j, are defined to be the majority spin and minority spin of 
the ferromagnet, respectively. 

The QD-electrode coupling is assumed to be sufficiently weak, thus, the sequential tun¬ 
neling is dominant and can be well described by the quantum master equation of reduced 
density matrix spanned by the eigenstates of the QD. The particle-number-resolved TCL 
quantum master equation for the reduced density matrix of the considered single QD is 
given by 


p'"' (i) 


dot,l 


-ic/"> - Y, 

aa 

+ Hx. 


( 4 ) 


For more details, see Methods section. Here, the complete basis {|0,0), |t;0), 14-, 0), 
is chosen to describe the electronic states of this single QD system, and the single QD system 
parameters are chosen as = ei = 1, U = 5, p = 0.9 and ksT = 0.04. 

Figure 2 shows the first four current cumulants as a function of the bias voltage for 
different ratios F^/F/j describing the left-right asymmetry of the QD-electrode coupling. 
We found that the non-Markovian effect has no influence on the current noise behaviors of 
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the single QD considered here, see Fig. 2. Scrutinizing Eq. Q, it is found that for the 
non-Markovian case the elements of the reduced density matrix are equivalent to that for 
the Markovian case because there are not the off-diagonal elements of the reduced density 
matrix. Thus, the equations of motion of the four elements of the reduced density matrix 
can be expressed as 


(0,0|P5 (t)|0,0) 

= - [rLt/i,+ (<^t) + ^RrfR,+ (^t)] (0)01 Ps (t) |o, 0) 

- [^LifL,+ + TRifR^+ (e^)] (0, 0| ps (t) |0, 0) 

+ [TL^h,- (6t) + (e^) (0, t| Ps (t) It, 0) 

+ [TlJl,- (6j + (ej (t, 0| ps (t) |0, t), (5) 


(0,t|p5 (t)|t,0) 

= [^LtfL,+ (et) + ^RrfR,+ (^t) (0, 0| Ps (t) |0, 0) 

- [^L-rfL- (et) + ^RtfR - (et)] (0, tl Ps (t) It, 0) 

- [^LifL,+ (et.^ - et) + ^RifR,+ (ett - e-^)] (0, tl Ps (t) It, 0) 

+ [r Lifi- (ett - et) + r RifR- (e^^ - e^) (t, 11 Ps (^) 11, t), (6) 

(t,0|ps (t)|0,t) 

= [r^t/T+ (et) + rRifR,+ (et) (0, 0| ps (t) |0, 0) 

- [^LifL- (et) + ^RifR- (et)] (t, 0| Ps (t) |0, t) 

- [rLt/L,+ (ett “ ^t) + ^RtfR,+ (ett “ <^t)] (t, 0| Ps {t) |0, t) 

+ (ett ~ ^t) + Ti^t/i?,- (ett “ ^^t) (t, tl Ps (t) |t,t), (7) 
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{iA\psit)n,i) 

= [^LifL,+ (e^ - et) + ^RlfR,+ “ '^t) (0>tl Ps (t) It, 0) 

+ [^LtfL,+ (et^ - e^.) + TR^fR^+ (e-^^ - e^) e“*^] (t, 0| ps (t) |0, t) 

- [^LtfL- (en ~ ^4^) + ^mfR - (^ti ~ ‘' 4 )] (t,tl Ps (t) It, t) 

- [^LifL- (et4 “ ^t) + ^RifR- (^44 “ ^^t)] (t,tl Ps (t) It, t) • (8) 

Here, fa,+ is the Fermi function of the electrode a, and fa- = 1 — fa,+ - The detailed 
procedure for calculation of the equation of motion of a reduced density matrix, see Methods 
section. Within the framework of the density matrix theory, the off-diagonal elements of the 
reduced density matrix characterize the qnantnm coherence of the considered QD system. 
Thus, the influence of the non-Markovian effect on the FCS may be associated with the 
quantnm coherence of the considered QD system. In order to confirm this conclusion, 
we take serially conpled and side-conpled donble QDs for illustration in the following two 
snbsection. 


Serially coupled double quantum dots with high quantum coherence 


We now consider two serially coupled double QDs weakly connected to two metallic 
electrodes, see Fig. 1(a). For the sake of simplicity, the spin degree of freedom has not been 
considered. The donble-QD is described by a spinless Hamiltonian 


H, 


dot, 2 




(9) 


where d\ {di) creates (annihilates) an electron with energy (which can be tnned by a gate 
voltage Vg) in ith QD. U is the interdot Coulomb repulsion between two electrons in the 
donble QD system, where we consider the intradot Coulomb interaction U —)■ cx), so that the 
donble-electron occupation in the same QD is prohibited. The last term of Hdot describes 
the hopping conpling between the two dots with J being the hopping parameter. The 
two metallic electrodes are modeled as non-interacting Fermi gases and the corresponding 
Hamiltonians are given by 


-^Leadsg E 

ak 


( 10 ) 
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where (flak) creates (annihilates) an electron with energy Sak and momentum k in a 
(a = L, R) electrode. The tunneling between the double QDs and the two electrodes is 
described by 

Ht ,2 = ^ (tLa\^^dl + + H.c. j . (11) 

ak 

For the case of the weak QD-electrode coupling, the particle-number-resolved TCL quantum 
master equation for the reduced density matrix of the considered serially double-QD system 
reads 


(() - (*) + p*”' (t) 

4piFp'”' (i) + p'”' (i) (t) d\ 

(«)4 

-4p‘”+‘’ (t) A'j^> + H.c. 


( 12 ) 


Here, we can diagonalize the serially coupled double QDs Hamiltonian Frdot ,2 in fhe basis 
represented by the electron occupation numbers in the QD-1 and QD-2 denoted respec¬ 
tively by Ni and Nj^, namely, {|0, 0), |1, 0), |0,1), |1,1)}, and obtain the corresponding four 
eigenstates of the considered serially coupled double QDs system [H] 

^dot,2 |0) = 0, |0) = |0,0), 

^dot,2 |1)^ = e± |1)^ , |1)^ = a± |1, 0) + |0,1), (13) 

.Hdot,2 |2) = 62 |2) = (ei -|- 62 + D) |2), |2) = |1,1), 


with 


e± 


(6i -F 6i) ± (6i - eif + 4J2 

2 


(14) 
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and 


a± 


b± 


tJ 

\j (e± — 

± (ezfc - ei) 

\l (e± — 


(15) 


Here, we focus on the regime (e+ — e_) k^T, where the hopping coupling between the two 
QDs strongly modihes the internal dynamics, and the off-diagonal elements of the reduced 
density matrix play an essential role in the electron tunneling processes [23l H91 - I5T] . In the 
following numerical calculations, thus, the parameters of the serially coupled double QDs 
system are chosen as ei = 62 = 1, J = 0.001, U = A and ksT = 0.05. 

When the coupling of the QD-2 with the right (drain) electrode is stronger than that of 
the QD-1 with the left (source) electrode, namely, < 1, we plot the hrst four current 

cumulants as a function of the bias voltage for different values of the QD-2-electrode coupling 
Fii; at Tl/Tr = 0.1 in Figs. 3(a)-3(d). We found that the non-Markovian effect has a very 
weak influence on the FCS. Interestingly, the high-order current cumulants the skewness 
and the kurtosis can still show the tiny differences, see Figs. 3(c) and 3(d). Whereas for 
the Ti/Tr > 1 case, the non-Markovian effect has a signihcant impact on the FCS, see Fig. 
4. Especially, for a relatively large value of the ratio Tl/Tr = 10 and the coupling of the 
QD-1 with the left electrode being stronger than the hoping coupling, namely, Tl/J > 1, 
the non-Markovian effect can induce a strong negative differential conductance (NDC) and 
super-Poissonian noise, see Figs. 4(e) and 4(f). In addition, in the case of Fl/F/j > 1 
and Tl/J > 1 , the transitions of the skewness and the kurtosis from positive (negative) to 
negative (positive) values are observed, see the dotted line in Fig. 4(c), the dotted and dash- 
dot-dotted lines in Fig. 4(d), and the dash-dot-dotted line in Fig. 4(h). It is well known 
that the skewness and the kurtosis (both its magnitude and sign) characterize, respectively, 
the asymmetry of and the peakedness of the distribution around the average transferred- 
electron number n during a time interval t, thus that provides further information for the 
counting statistics beyond the shot noise. 

To discuss the underlying mechanisms of the current noise clearly, for the system param¬ 
eters considered here, the two singly-occupied eigenstates and eigenvalues can be expressed 
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as 


|l)^ = T#|l,0) + f |0,1), 


2 I-*-’ ' 2 

€+ = €_ = e. 


(16) 


Here we have utilized the equations ei = €2 = e and e 3> J. In this situation, the equations 
of motion of the six elements of the reduced density matrix are given by 


(0,0|ps(«)|0,0) 

= - [ri/i,+ (e) + r„/„.+ (e)] (0,0| PS (i) |0,0) 

+ i [Fsh.- (£) + rgfn,. (e) e‘^] (11+ ps («) |1)+ 

- i [ri/p_ (£) - r«/H,_ (€) eq (1|+ PS («) |1)- 

- I [U/i,- (£) - rfl/R,_ (€) e'x] (ir PS («) |1)+ 

+ 2 (e) + Fr/h,-(€) e'^] (1| ps{t)\l) , 


( 17 ) 


(11*4”’ (*)|1)* 

= ^ [ri/i,+ (e) + Fr/rr (e) e-«] (0,0| ps (t) |0,0) 


1 


Tq, [fa,+ {e + U) + fa - (e)] (1|^ Ps (t) |1)’ 


a=L,R 

± (iti =F xFr) (1|+ PS (t) |1)- T (its T rFg) (1|+ ps (t) |1)- 


ir. 


ir 


± ttFs) (ir PS (t) ii)^ ±~ ± (ir PS (t) ii)- 

+ ^ [FsA- (e + f/) + (e + [/) (1,11 ps (t) |1,1), 


( 18 ) 


9 



{ii*w(i)iir 

= -l [rL/L,+ (e) - r„/„.+ (e) e-'*] (0,0| ps (i) |0.0) 

± T rfi) (ir pSt’ («) ii)* =F (»H =F (ir PS («) ii>-" 

- ^ E [/o,+ (f + c^) + A- W1 (1|*PS(«) |1>^ T 2a(i|±ps (i) |1)^ 

a=L^R 

T (*4i ± ttF,) (ir pf («) |1)" ± (**« ± (ir Ps («) |1>" 

+ ^ [ri/i,- (6 + C/) - r„/H,_ (e + t/) e«] (1,1| PS («) |1,1>, (19) 

{l,l|ps(«)|l,l) 

= I [ri/i,+ (e + C/) + rs/s,+ (6 + C/) e-‘*] {1|+ PS («) |1>+ 

+ I [ri/i,+ (e + C/) - rs/H,+ (e + U) e-«] (l|+ps (i) 11)’ 

+ ^ [rL/i,+ (e + C/) - r„/«,+ (e + U) e-*] (ip Ps (i) |1)+ 

+ 2 [1's/l,+(^ + fl) + rR/fi,+(^ + fl)s {1| Ps(l)|l) 

- iri/i,_(€ + (7) + rR/R,_(£ + (7)](i,i|ps(«)|i.i). ( 20 ) 

where $q, = (e + t/) — 0^ (e), 4>a (e) = Re^ | + ^ 2 Td!^T digamma function) 

and Fa = fa,+ i^ + U) — fa - (e). Compared with the Markovian case, it is obvious that 
the non-Markovian effect manifests itself through the off-diagonal elements of the reduced 
density matrix, namely, the quantum coherence of the considered QDs system. In Fig. 
5(a), we plot the functions — O.ld’j? (!"/? = O.IFi), (r_R = F^,) and 0.1<h2, — <Fr 

(Fi = O.lF/j) as a function of bias voltage. It is clearly evident that the values of the 
functions — 0.1<h/{ and <hz, — <h_R show significant variations with increasing bias voltage, 
especially in the vicinity of the bias voltages H = 2 and 14 = 10 because the new transport 
channels begin to participate in quantum transport; while 0.1$^ — has a gentle variation. 
Consequently, the non-Markovian effects in the F^/F/j > 1 case have a remarkable impact 
on the FCS, see Fig. 4. Moreover, for F^/F^ = 10 case, the non-Markovian effect has a more 
signihcant on the FCS than the F^/Fr = 1 case, which originates from the QD-2-electrode 
coupling Fk is weaker than the hoping coupling J, where the electron tunneling from QD-1 
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can not tunnel out QD-2 very quickly and still influence the internal dynamics. 

In order to illustrate whether the non-Markovian effect has a weak influence on the PCS 
in a relatively small quantum-coherent QD system, we consider the regime (e+ — e_) 3> ksT 
(J = 1), where the off-diagonal elements of the reduced density matrix have little influence 
on the electron tunneling processes. We And that for the J = 1 case the diagonal elements 
of the reduced density matrix play a major role in the electron tunneling processes, and the 
non-Markovian effect in this case indeed has little impact on the PCS, see Pigs. 3(e)-3(h). 
Consequently, the influence of the non-Markovian effect on the PCS depends on the quantum 
coherence of the considered QD system. To prove whether this conclusion is universal or 
not, we take side-coupled double QDs for further illustration in the following subsection. 


Side-coupled double quantum dots with high quantum coherence 

We consider here a side-coupled double QDs system. In this case, the QD-1 is only weakly 
coupled to the two electrodes, see Pig. 1(b). The QD-electrode tunneling is thus described 
by 

Ht, 3 = (^ak^akdi + • ( 21 ) 

In the case of the QD-electrode weak coupling, the particle-number-resolved TCP quan¬ 
tum master equation for the side-coupled double QDs can be expressed as 


L 

+p<”> (() 


-1 w 


-4p*”+'’(«).4S,+’ + H.c. . 


( 22 ) 


Here, the eigenstates and eigenvalues of the side-coupled double QDs system are the same 
as the serially coupled double QDs system. In the following numerical calculations, the 
parameters of the side-coupled QDs system are chosen as ei = 62 = 1, J = 0.001, U = 5 and 
keT = 0 . 1 . 

Por the present side-coupled QDs system with high quantum coherence, we And that for 
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Tl/F/j > 1 case the non-Markovian effect has a more remarkable impact on the FCS than 
that in the serially coupled double QDs system, but the NDC does not appear, see Figs. 4 
and 6. For instance, in the case of F^/J > 1 and Tl/Tr = 1, the non-Markovian effect can 
further enhance the super-Poissonian shot noise, see the dotted and dash-dot-dotted lines 
in Fig. 6(b); and the transitions of the skewness and the kurtosis from a relatively small 
positive to a large negative values take place, especially for a relatively large value F^^/J 
the kurtosis can be further decreased to a very large negative value, see the dotted and 
dash-dot-dotted lines in Figs. 6(c) and 6(d). While for the F^/J > 1 and F^/F^j = 10 case 
the non-Markovian effect can enhance the shot noise to a super-Poissonian value, see the 
dotted and dash-dot-dotted lines in Fig. 6(f), and the transition of the kurtosis from small 
positive to large negative values only takes place, see the dotted and dash-dot-dotted lines 
in Fig. 6(h). For the system parameters considered here, namely, in the limit of e 3> J, the 
equations of motion of the six elements of the reduced density matrix read 

(O.OIpsW |0,0> 

= - |ri/i,+ (t) + rK/H,+ (t)] (0, o| PS it) |o, 0) 

+ ) [rph.- (€) + r„/K,_ (£)e«] (ir |i)+ 

- 2 (e) + Th/h,- (€)eq (1|^ Ps (^) 11) 

- I [Ft/i,-(€) + rK/„,_ (e)e"‘] (ir ps(«) |1>+ 

+ ) [ri/p- (€) + r„/s.- (e) e-x] (ir PS (() |1)- , (23) 
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= ^ [rLfL,+ (e) + Tnfn,+ (e) (0,0| ps (t) |0, 0) 

- ^ E (^ + U) + /«,- (e)] (1|^ PS it) |1)^ 

a=L,R 

a=L,R 

=F^ E ^l*4«±>r-K,]{l|>s(«)|l>+ 


+ 


1 

2 


TLfL,-{e + U) + hRfR,_ {e + U)e^^ 


( 1 , 1 |P 5 (t)|l,l), 


(24) 


(l|E 5 (t)|ir 

= -l [rLfL,+ (e) + r«/ij,+ (e) (0, 0| ps (t) |0, 0) 

a=L,R 

- ^ E r« [/«.+ (f + C') + /«,- W1 (i|^ PS (t) |i>^ T 2a {i|±ps (i) |i)T 

a=L,R 


T 


1 

2 


a=L,R 


+ I [Fl/l,- (6 + C/) + r„/fl,_ (6 + tJ) e'*] (1,11 PS (i) |1.1), 


(25) 


(l,l|ps(«)|l.l) 

= I [ri/i,+ (6 + c/) + rH/„,+ (e + u) €-'<] { 1 |+ PS («) | 1 >+ 

+ i [ri/L+ (e + C/) + rH/R,+ (e + U) e-‘x] (11+ ps («) ll)' 

+ ^ [ri/i,+ (e + C/) + rH/s,+ (e + U) e-«] (Ip ps («) |1>-" 

+ ) [rL/i,+ (e + C/) + rB/R+ (e + C/) e-‘x] (1 p ps («) 1 1>' 

- [ri/i,_ (6 + t/) + Fr/h,. (e + t/)] (1,1| PS («) |1, 1). (26) 

From the above four equations, we find that these characteristics also originate from the 
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quantum coherence of the side-coupled double QDs, and can also be understood in terms of 
the functions and which have considerable variations in the vicinity of 

the bias voltages Vb = 2 and H = 12 because the new transport channels begin to enter the 
bias voltage window, see the solid and dashed lines in Fig. 5(b). As for the Vl/Vr < 1 case 
the non-Markovian effect has a slightly influence on the FCS because the function 
has a gentle variation with increasing the bias voltage, see the dotted line in Fig. 5(b), which 
is the same as the serially coupled double QDs system, see Figs. 3(a)-3(d) and 7. 

In addition, it should be pointed out that for Ti/Tr = 1 the non-Markovian effect has 
a stronger impact on the FCS than that for Tr/Tr > 1 case, which is contrary to the case 
of the serially coupled double QDs system. For the the side-coupled double QDs system, 
the quantum coherence originates from the quantum interference between the direct electron 
tunneling process, namely, the conduction-electron tunneling into the QD-1 and then directly 
tunneling out of the QD-1 onto the drain electrode, and the indirect tunneling process, 
namely, the conduction-electron from the source electrode first tunneling from the QD-1 to 
the QD-2, then tunneling back into the QD-1, and at last tunneling out of the QD-1. Thus, 
the fast direct tunneling process in the F 2 , = lOF/? case can be suppressed compared with 
the Fi, = F/j case, which leads to the non-Markovian effect has a relatively strong impact 
on the FCS in the Vl/Vr = 1 case. 

Discussion 

We have developed a non-Markovian FCS formalism based on the exact TCL master 
equation, and studied the influence of the interplay between the quantum coherence and 
non-Markovian effect on the long-time limit of the FCS in three QD systems, namely, single 
QD, serially coupled double QDs and side-coupled double QDs. It is demonstrated that the 
non-Markovian effect manifests itself through the quantum coherence of the considered QD 
molecule system, and especially has a signihcant impact on the FCS in the high quantum- 
coherent QD molecule system, which depends on the coupling of the considered QD molecule 
system with the source and drain electrodes. For the single QD system without quantum 
coherence, the non-Markovian effect has no influence on the current noise properties; whereas 
for the serially coupled and side-coupled double QDs systems with high quantum coherence, 
that has a remarkable impact on the FCS when the coupling of the considered QD molecule 
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with the incident electrode is equal to or stronger than that with the outgoing electrode. 
For instance, for the high quantum-coherent serially coupled double QDs system, the non- 
Markovian effect can induce a strong NDC and change the shot noise from the sub-Poissonian 
to super-Poissonian distribution in the case of Pl/Pj? ^ 1 and P^ > J; while for the 
high quantum-coherent side-coupled double QDs system, that can remarkably enhance the 
super-Poissonian noise or the sub-Poissonian noise for the Tl/Tr > 1 case. Moreover, the 
non-Markovian effect can also lead to the occurrences of the skewness and kurtosis from 
small positive to large negative values. These results indicated that the influence of the 
non-Markovian effect on the long-time limit of the PCS should be considered in a highly 
quantum-coherent single-molecule system. 

Methods 

Particle-number-resolved time-convolutionless quantum master equation 

We consider a general transport setup consisting of a single-level QD molecule weakly 
coupled to the two electrodes, see Fig. 1, which is described by the following Hamiltonian 

H = hfelectrodes T hfdot T .^hyb- (2'^) 

Here, the hrst term hfeiectrodes = Xla k a stands for the Hamiltonians of the two 

electrodes, with Sak being the energy dispersion, and ttaka (^Lfco-) annihilation (creation) 
operators in the a electrode. The second term iJ^ot = Us {dl., d^), which may contain vibra¬ 
tional or spin degrees of freedom and different types of many-body interaction, represents the 
QD molecule Hamiltonian, where (d^) is the creation (annihilation) operator of electrons 
in a quantum state denoted by /i. The third term ddhyb = 

describes the tunneling coupling between the QD molecule and the two electrodes, which is 
assumed to be a sum of bilinear terms that each create an electron in the QD molecule and 
annihilate one in the electrodes or vice versa. 

The QD-electrode coupling is assumed to be sufficiently weak, so that ddhyb can be treated 
perturbatively. In the interaction representation, the equation of motion for the total density 
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matrix reads 


(«) = -“ [-ffhjb (<). p' (*)]=£(<) p' (i). 


( 28 ) 


with 

^hyb it) = -Yl 4 

a,fi 

where (t) = and d^ (t) = In order to 

derive an exact eqnation of motion for the rednced density matrix ps of the QD molecule 
system, it is convenient to define a super-operator V according to 


Vp = trs [p] ® = p5 ® p^, 


(29) 


with pb being some fixed state of the electron electrode. Accordingly, a complementary 
super-operator Q reads 


Qp = p- Vp. (30) 

For a factorizing initial condition p(fo) = Ps (to) ® Pb, Vpit^) = pito), and Qp(fo) = 0. 
Using the TCL projection operator method 1^2], one can obtain the second-order TCL 
master equation 

= dhVCit)C{h)Vp{t), (31) 


The Eq. (31) is the starting point of deriving the particle-number-resolved quantum master 


equation. Using Eqs. (28) and (29), after some algebraic calculations we can rewrite Eq. 


(31) as 


d 

si"-" (‘) 

- / dtitiB [pi,s {t) 0 PbIL (ti) dy (fi) dl [t) (t)] 

afiv d-oc 

~ K (^) (^) PPS {t) ® Pb] 

afiu d-oo 

+ ^ f dtitiB [fil, (t) d^ (t) Pi^s (t) ® pBdl (tl) fau (U)] 

+ ^ [ dtitiB [dl {t) fai, (t) PI,S {t) 0 PbJL (ti) du (U)] + H.C. 


(32) 


afiu 


16 






In order to fully describe the electron transport problem, we should record the number of 
electrons arriving at the drain electrode, which emitted from the source electrode and passing 
through the QD molecule. We follow Li and co-authors [531IM] and introduce the Hilbert 
subspace {n = 1,2,...) corresponding to n electrons arriving at the drain electrode, 
which is spanned by the product of all many-particle states of the two isolated electrodes, 
and formally denoted as = span| 0 Then, the entire Hilbert space of 

the two electrodes can be expressed as B = With this classihcation of the electrode 


states, the average over states in the entire Hilbert space B in Eq. (32) should be replaced 


with the states in the subspace and leading to a conditional TCL master equation 

(t) 

dtPi,s 

= -^ [ [pi^s {t) ® PsfL ih) (H) dl (t) (t)] 

a^iu d-oo 

[ dtitiBin) [dl {t) (t) fl^ (ti) (ti) pps {t) 0 Pb] 

a^iu d-oo 

+ '^ [ dtitr^M [fl^ (ti) d^, {ti) pps {t) 0 psdl {t) /„/, {t)] 

aiiu d-oo 

+ ^ [ dtitiBCu) [dl (t) faf, (t) pps (t) ® PsfL (ti) du (H)] + H.C.. (33) 


afii/ 


To proceed, two physical considerations are further implemented, (i) Instead of the 
conventional Born approximation for the entire density matrix pr {t) — p (t) 0 ps, the 
ansatz p^ (t) ~ p^"^ (t) 0 p^^ is proposed, where p^^ being the electrode density operator 
associated with n electrons arriving at the drain electrode. With this ansatz for the entire 
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density operator, tracing over the snbspace the Eq. (33) can be rewritten as 


— it) 

dtPi,s 

“ / dtitiBin) [fl^ (ti) {t) pb] pPs (t) du ih) dl {t) 

afiu 

“ / dtitiBin) [fa^ {t) fl^ (ti) Pb] dl (t) d^, {h) pf]. {t) 

afiu d-oo 

+ 5Z / dtitiBCn) pB d^{ti) p^^l{t)dl{t) 

fiu d-oo L -I 

+ 5Z / dtitTsi„) Pb d^{B)p^^~^\t)dl{t) 

^lU d-oo 

+ 5Z / dtitiBin) fl^ih) fL,,{t) Pb dl{t) p^^l{t)d^,{ti) 

^lU d-oo L -I 

+ 5^/ d^itr^w flAti) Pb dl{t) du {h) + H.c.. 


flip 


(34) 


Here we have used the orthogonality between the states in different subspaces, (ii) The 
extra electrons arriving at the drain electrode will ffow back into the source electrode via the 
external closed transport circuit. Moreover, the rapid relaxation processes in the electrodes 
will bring the electrodes to the local thermal equilibrium states quickly, which are determined 


by the chemical potentials. Consequently, after the procedure done in Eq. (34), the electrode 


density matrices and should be replaced by p^\ In the Schrodinger representation. 
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the Eq. (34) can be expressed as 


4„i") {f. 


= -i [-ffs, 4”’ (<) 

aiJ.u 

aflu 

flU 

+ Y f *i4;l, (« - «i) (i)< 

flU 

+ ^ / (ti - t) (*) + H.c. 


(36) 


/iZ/ 


where the correlation function are dehned as 


Ci^Ju trR [fl^ (t) fau ih) Pb] = {fif, (t) fau (ti)> , 

Ci~l it-ti)= trjj (t) fl^ (ti) Pb] = {fatp (t) fl^ (ti)) . (36) 


Introducing the following super-operators 

u 

(<) = E /' “'‘■C'U - ‘i) (37) 
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then, the Eq. (35) can be rewritten as a compact form 


(+) 


* [Hs, (t) 

Y1 ^ (^) pP (^) 

(*) ^>s ’ (*) <*1. “ (<) Ps”^‘’ (<) <*!. 

(‘) (‘) - 4p!r*‘’ (‘) pi£ (‘) + . 


(38) 

where A^p {t) = J2<x (t). The above equation is the starting point of the non-Markovian 
FCS calculation. 


Non-Markovian full counting statistics 


In this subsection, we outline the procedure to calculate the non-Markovian FCS based 


on Eq. (38). The FCS can be obtained from the cumulant generating function (CGF) F (x) 
which related to the probability distribution P{n,t) by [Ml [55] {n,t) 

where y is the counting held. The CGF F (y) connects with the particle-number- 
resolved density matrix (t) by dehning S (y,t) = (t) e*"^. Evidently, we have 

g--^(x) =Tr[S' (y, t)], where the trace is over the eigenstates of the QD molecule system. Since 


Eq. (38) has the following form pG) = Ap^'^'> + + Dip^'^ then, S (y,t) satishes 


S = AS + e *^(7S' + F^DS = T^S*, where S' is a column matrix, and A, C and D are three 
square matrices. The specihc form of can be obtained by performing a discrete Fourier 


transformation to the matrix element of Eq. (38). In the low frequency limit, the counting 


time, namely, the time of measurement is much longer than the time of tunneling through 
the QD molecule system. In this case, F (y) is given byjSl llOl SSI jSSHSTj F (y) = — Ai (y) t, 
where Ai (y) is the eigenvalue of which goes to zero for y —)• 0. According to the dehnition 
of the cumulants one can express Ai (y) as Ai (y) = YlT=i order cumulants 

can be calculated by the Rayleigh-Schrodinger perturbation theory in the counting param¬ 
eter y. In order to calculate the hrst four current cumulants we expand to four order in 
X 


Lx — Lq + Lix + XyL2X + XfLsX + TiLiX + 


2 ! 


3!' 


4!' 


(39) 
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and define the two projectors [IHl SSI ESI EH] P = = |0)) ((O | and Q = Q‘^ = 1 — P, 

obeying the relations PLq = LqP = 0 and QLq = LqQ = Lq. Here, |0)) is the right 
eigenvector of Lq, i.e., Lq |0)) = 0, and ((o| = 1 is the corresponding left eigenvector. In 
view of Lq being regnlar, we also introdnce the psendoinverse according to R = 
which is well-dehned dne to the inversion being performed only in the snbspace spanned by 
Q. After a carefnl calcnlation, Ai (y) is given by 

^1 (x) = ((0|Li|0))x 

+ ^_[{{0\L,\0))-2{{0\L,RL,\0))] 

+ ^ [((0 I Ls 10)) - 3 ((0 I {L^RL, + L.RL^) \ 0)) 

-6 ((0 I L,R {RLiP - LiR) h \ 0))] y=*+ 

+ 1_[{{0\L,\0))-6{{0\L,RL,\0)) 

-4{{0\{LsRL, + L,RLs)\0)) 

- 12 ((0 1 L 2 R {RL^P - LiR) Li 1 0)) 

- 12 ((0 I LiR {RL 2 P - L 2 R) Li 10)) 

- 12 ((0 I LiR {RLiP - LiR) L 2 10)) 

- 24 ((0 I LiR (i^^LiPLiP - RLiPLiR - LiR^LiP 

-RLiRLiP + LiPLiP) Li 10))] / + • ■ ■ . (40) 


From Eq. (40) we can identify the hrst fonr cnrrent cnmulants: 


Ci/f = ((0|Li|0))A, 


C2/t = [((0 I L2 10)) - 2 ((0 I L,RL, 10))] /^^ 


(41) 

(42) 


C,/t = [((0 I L3 10)) - 3 ((0 I (L2PL1 + L1PL2) 10)) 

-6 ((0 I LiR (PLiP - LiP) Li 10))] /i\ (43) 
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C^jt = [((0 I L 4 | 0 )) - 6 ((0 I URU 1 0 )) 

-4((0|(L3i?Li + Lii?L 3 )| 0 )) 

- 12 ((0 I UR {RUP - UR) U 1 0 )) 

- 12 {{0\ UR{RUP - UR) U | 0 )) 

- 12 (( 0 1 UR {PUP - UR) U 1 0 )) 

- 24 (( 0 1 UR {R^UPUP - RUPUR - UR^UP 

-RURLiP + URUR)U\0))]/i\ (44) 

Here, it is important to emphasize that the first four cumulants Ck are directly related to the 
transport characteristics. For example, the hrst-order cumulant (the peak position of the 
distribution of transferred-electron number) Ci = n gives the average current (/) = eCi/t. 
The zero-frequency shot noise is related to the second-order cumulant (the peak-width of 
the distribution) S = 2 e^C 2 /t = 2e^ {u — /t. The third-order cumulant C 3 = (n — n)^ 

and four-order cumulant C4 = (n — n) —3{n — n) characterize, respectively, the skewness 
and kurtosis of the distribution. Here, (■ ■ ■) = Xln i''^ general, the shot noise, 

skewness and kurtosis are represented by the Fano factor F 2 = C 2 /C 1 , F 3 = C^/Ci and 
F 4 = C 4 /C 1 , respectively. 
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(a) 


Left electrode ^ ^ Right electrode 




FIG. 1: Schematic of the two single-level QD molecules weakly coupled to two electrodes, (a) 
serially coupled double QDs, (b) side-coupled double QDs. Here, the two QD molecules possess high 
quantum coherence in the case of A ^ ksT (A being the singly-occupied eigenenergy separation, 
ks the Boltzmann constant, T the temperature of the QDs system). The hopping coupling between 
the two QDs, and the strength of coupling between the QDs system and the electrode a, 
characterized by J and Fq, respectively. 


are 









FIG. 2; The average current ((/)), shot noise (C' 2 /C'i), skewness {C^/Ci) and kurtosis {Ca^/Ci) 
versus bias voltage for the Morkovian and the non-Markovian case at different coupling of the 
single QD with two ferromagnetic electrodes, respectively. Here, Ck is the zero-frequency fc-order 
cumulant of current fluctuations. The non-Markovian effect has no influence on the first four 
current cumulants of the considered single QD. The single QD system parameters: e^i- = = 1, 

U = 5, p = 0.9 and k^T = 0.04. 
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FIG. 3; (a)-(d) The average current ((/)), shot noise {C 2 /C 1 ), skewness (Ca/Ci) and kurtosis 
(C' 4 /C'i) versus bias voltage for the Morkovian and the non-Markovian case at different values of the 
QD-2-electrode coupling F^j with F l/^r = 0.1. Here, Ck is the zero-frequency /c-order cumulant of 
current fluctuations. The non-Markovian effect in the T^/F/j = 0.1 case has a weak influence on the 
the first four current cumulants. The serially coupled double QDs system parameters: ei = 62 = 1, 
J = 0.001, 1/ = 4 and k^T = 0.05. (e)-(h) The average current ((/)), shot noise {C 2 /C 1 ), skewness 
{C 3 /C 1 ) and kurtosis (C' 4 /C'i) versus bias voltage for different coupling of the serially coupled 
double QDs system with two metallic electrodes. Here three cases are considered, namely, (1) the 
Markovian and the diagonal elements of the reduced density matrix, (2) the Markovian and the 
off-diagonal elements of the reduced density matrix, (3) the non-Markovian and the off-diagonal 
elements of the reduced density matrix. The non-Markovian effect has a very weak influence on 
the first four current cumulants in the serially coupled double QD system with a relatively small 
quantum coherence. The serially coupled double QDs system parameters: ei = 62 = 1, J = 1, 
U = A and ksT = 0.05. 
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FIG. 4: The average current ((/)), shot noise (C' 2 /C'i), skewness (Ca/Ci) and kurtosis {Ca^/Ci) 
versus bias voltage for the Morkovian and the non-Markovian case at different values of the QD- 
2-electrode coupling Fij. (a)-(d) for F/^/Fij = 1, (e)-(h) for F^/F/^ = 10. Here, Ck is the zero- 
frequency fc-order cumulant of current fluctuations. The non-Markovian effect in the Tx/F^ > 1 
case has a significant impact on the first four cumulants of transport current. The serially coupled 
double QDs system parameters: ei = £2 = 1, J = 0.001, U = A and ksT = 0.05. 
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FIG. 5: (a) The functions (r/j = 0.iri), (ri{ = r 2 ,) and 

(Fi = O.lF/j) as a function of bias voltage with U = A and ksT = 0.05. (b) The functions 

‘hi + 0.1<I>/j (Fij = O.IF l), (Fij = F^,) and 0.1$^ + ^r (F^, = O.IFr) as a function of bias 


voltage with U = 5 and ksT = 0.1. Here, <I>c 


= c/)a {e + U) - (t)a (e), (e) = Re T \ + 


and T is the digamma function. The variation of the value of the above mentioned function 
is responsible for whether the non-Markovian effect has a remarkable influence on the first four 
cumulants of transport current. 
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FIG. 6: The average current ((/)), shot noise (C' 2 /C'i), skewness (Ca/Ci) and kurtosis {Ca^/Ci) 
versus bias voltage for the Morkovian and the non-Markovian case at different values of the QD- 
1-electrode coupling Fij. (a)-(d) for Ti/Tn = 1, (e)-(h) for F^/F/^ = 10. Here, Ck is the zero- 
frequency fc-order cumulant of current fluctuations. The non-Markovian effect in the Tx/Fij > 1 
case has a more remarkable impact on the first four cumulants of transport current than that in 
the serially coupled double QDs system, but the NDC does not appear. The side-coupled double 
QDs system parameters: ei = 62 = 1, J = 0.001, U = 5 and ksT = 0.1. 
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FIG. 7: The average current ((/)), shot noise (C' 2 /C'i), skewness (Ca/Ci) and kurtosis {Ca^/Ci) 
versus bias voltage for the Morkovian and the non-Markovian case at different values of the QD- 
1-electrode coupling F^j with Ti/Tr = 0.1. Here, Ck is the zero-frequency A:-order cumulant of 
current fluctuations. The non-Markovian effect in the Fi/Fjj = 0.1 case has a slightly influence on 
the the first four current cumulants. The other system parameters are the same as in Fig. 6. 
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